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We consider the nonequilibrium evolution of an O(N)— symmetric scalar quantum 
field theory using a systematic two-particle irreducible 1/N— expansion to next-to- 



\J^) ' leading order, which includes scattering and memory effects. The corresponding 



"full Kadanoff-Baym equations" are solved numerically without further approx- 
imations. This allows one to obtain a controlled nonperturbative description of 
1 far-from-equilibrium dynamics and the late-time approach to quantum thermal 

^— ^ ' equilibrium. Employing in addition a first-order gradient expansion for the Wigner 

transformed correlators we derive kinetic equations. In contrast to standard de- 
scriptions based on loop expansions, our equations remain valid for nonperturba- 
tively large fluctuations. As an application, we discuss the fluctuation dominated 
regime following parametric resonance in quantum field theory. 



In recent years we have witnessed an enormous increase of interest in 
the dynamics of quantum fields out of equilibrium. Strong motivation in 
elementary particle physics comes, in particular, from current and upcom- 
ing relativistic heavy-ion collision experiments, phase transitions in the 
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early universe or astrophysical applications. Here the involved nonequi- 
librium dynamics is often characterized by large corrections from quantum- 
statistical fluctuations that are not accessible in a weak coupling or loop 
expansion. A paradigm for such a situation is provided by the phenomenon 
of parametric resonance, which represents an important building block for 
our understanding of the (pre) heating of the early universe after a pe- 
riod of inflation. 1 In this context the resonant amplification of fluctuations 
leads to explosive particle production with a transition from a classical to 
a fluctuation dominated regime, characterized by nonperturbatively large 
occupation numbers inversely proportional to the coupling. 

Until recently 2 , classical field theory studies on the lattice have been 
the only quantitative approach available 3 . These are expected to be valid 
for not too late times, before the approach to quantum thermal equilibrium 
sets in. Calculations in quantum field theory had been limited to mean-field 
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type approximations (leading-order in large-TV, or Hartree) which neglect 
scatterings. 4 These approximations are known to fail to describe late-time 
thermalization and, even at early times, do not give a valid description of 
the entire amplification regime as pointed out in Refs. 2,5 . 

In Ref. 2 the first study of parametric resonance in quantum field theory 
from a next-to-leading order (NLO) calculation in a systematic two-particle 
irreducible (2PI) 1/iV-expansion 6 ' 7 was presented, which includes scatter- 
ing and memory effects. The classical resonant amplification at early times 
is found to be followed by a collective amplification regime with explo- 
sive particle production in a broad momentum range. As a consequence, 
one observes rapid prcthcrmalization with a particle number distribution 
monotonous in momentum. In particular, in this regime there are leading 
contributions from all 2PI loop orders and standard weak coupling or loop 
expansions break down. For its description it is crucial to employ a non- 
perturbative approximation as provided by the 1 /TV-expansion at NLO. 2 

In this note, after reviewing the nonperturbative physics involved in the 
phenomenon of parametric resonance, we present suitable kinetic equations 
derived 8 from the 2PI 1 /TV-expansion at NLO 6 . For similar approximation 
schemes see also Ref. 9 . Standard descriptions are typically based on a loop 
expansion of the 2PI effective action. In particular the classical Boltzmann 
equation can be obtained starting from a three-loop approximation. 10,11,12 
The advantage of our kinetic equations is that their applicability is not 
limited by the question of the validity of a loop expansion. For sufficiently 
large N, 13 it is only restricted by the applicability of a gradient expansion 
employed in the derivation of the kinetic equations. Apart from the ex- 
ample of the fluctuation dominated regime following parametric resonance, 
as a further important application these equations are expected to be valid 
near second-order phase transitions. Here the correlation length diverges at 
the transition in the static limit. The possibility of enhanced fluctuations 
near such a critical point for the nonequilibrium dynamics has received 
much attention recently in the context of relativistic heavy-ion collisions. 14 
A detailed discussion will be presented elsewhere. 8 

We consider a real scalar A-component quantum field ip a (a= 1, . . . , A) 
with A/(4!A r ) (ip a (p a ) 2 interaction, where summation over repeated indices 
is implied. All correlation functions of the quantum theory can be ob- 
tained from the 2PI generating functional for Green's functions T[(f), G], 
parametrized by the field expectation value <p a {x) = (<p a (x)) an d the con- 
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Figure 1. Total energy Etot (solid line) and classical field energy E c \ (dotted line) as a 
function of time for A = and = 4 in 3 + 1 dimensions, from Ref. 2 . The dashed 

line represents the fluctuation part, showing a transition between a classical field and a 
fluctuation dominated regime. For t > tnonpert there are leading contributions from all 
loop orders and the employed nonperturbative approximation becomes crucial. 



nected propagator G ab (x,y): 
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T[<f>, G] = S[tf>] + ^TrlnG- 1 + '-Tr G^(4>) G + T 2 [<j>, G] . 

Here iG$ ab (x,y; (f>) = 5 2 S[4>]/54>a{x)8<j) b {y), where S is the classical action 
with = — f </> a (n x + m 2 )(j) a /2 as the free part. We use the notation 
/ = J c dx° J dx with C denoting a closed time path along the real axis. 16 
At NLO in the 1/iV-expansion of the 2PI effective action one has 6,7 

r 2 [^,G] = - A J G aa (x,x)G bb (x,x) + ^Trln[B(G)] 



+ ^ J I{x,y,G)<j) a (x)G ab (x,y)<f> b (y) , (1) 

B(a;, y; G) = 8 c (x - y) + i-^G ab (x, y)G ab (x 1 y) , (2) 
l(x,y;G) = -^G ab (x,y)G ab (x,y) 



-i-j^f J r ( x ' z ; G)G ab (z, y)G ab {z, y) . 



(3) 



Here the function I resums an infinite number of "chain graphs" . 6 The 
equations of motion for (j) a and G ab are given by 15 

S_mG] =Q ST[^G] _ Q 

5<t> a (x) ' SG ab (x,y) 

The resulting evolution equations have been derived in detail in Refs. 6 ' 7 . 
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Figure 2. Effective particle number density for the transverse modes as a function of 
time for various momenta < p < 5po an d same parameters as in Fig. 1. At times 
t < tnonport there is exponential amplification, which affects a broad momentum range 
for t > t,£_ decay Taken from Ref. 2 . 



They have been solved numerically in Ref. 2 on a lattice along the lines of 
Ref. 6 . We decompose the full two-point function using 

G ab (x, y) = G> (x, y) 6 C (x° - y°) + G< (x, y) 6 C (y° - x°) (5) 

such that the spectral function p and the symmetric propagator F are given 
by 6 < 17 

Pab{x,y) = -2lm(G> b {x,y)^j and F ab (x, y) = Re(G> b (x, y) ) (6) 

For the discussion of parametric resonance in Ref. 2 a system is considered, 
which is initially in a pure quantum state . Employing spatially homoge- 
neous fields 4> a (t = t ) ~ (67V/A) 1 / 2 , the initial propagator is taken to be 
diagonal with F a b — diag{F||,FL, . . . ,F±} and equivalently for p a \,. The 
effective transverse particle number density 
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n±(p) + 77= F ± (t,t';p)d t d t ,F±(t,t';p) (7) 



1/2 



t=t' 



and corresponding longitudinal particle number density ray are defined as 
in Refs. 6,17 . Initially n± = nn =0 such that the initial conditions are 
characterized by zero particle number and large field amplitudes for small 
coupling A. Correspondingly, the conserved total energy E to t is initially well 
approximated by the classical field contribution, i.e. E to t — E c \, as shown in 
Fig. 1. When the system evolves in time, more and more energy is converted 
into fluctuations. A detailed discussion of the various characteristic regimes 
is given in Ref. 2 . For times t > i n0 npcrt this leads to nonperturbatively large 
particle number densities inversely proportional to the coupling. Neglecting 
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the difference between transverse and longitudinal modes at sufficiently late 
times, with F a b ~ F S a b and p a b ~ p8 a b we have 2 



(8) 



The nonperturbatively large enhancement of the statistical propagator F 
has the important consequence that any approximation scheme based on a 
loop expansion of the 2PI generating functional breaks down. We emphasize 
that the 2PI 1/iV-expansion at NLO remains valid as pointed out in Ref. 2 . 
To discuss this in more detail, we consider here the simplified evolution 
equations for F and p with cf) — 0. For the general case with a nonvanishing 
field expectation value see Ref. 2 . One finds 6 ' 17 



(9) 



[D x + M 2 {x)} F(x, v)=-f dz°j dz X p (x, z)F(z, y) 
+ / dz°/dz E F (x,z)p(z,y), 

Jt J 

p x + m 2 (x)\ P {x, y) = -J y0 dz °J dz s p( x ' z )p( z ' v)- 

At NLO in the 1 /iV-expansion the effective mass term M 2 (x) is given by 



MHx) = m + X- 



6N 



-F(x,x) 



and the self-energies are 6 

T, F (x,y) F{x,y)I F {x,y) - ^p(x,y)I p (x,y) 

Z p (x,y) p(x,y)I F (x,y) + F(x,y)I p (x,y) . 

Here the functions I F and I p contain the resummation: 

A A f x ° f 

I F (x,y) = --Il F (x,y) + - J dz" J dz I p (x,z)U F (z,y) 



(10) 

(11) 
(12) 



> to 
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dz I F (x,z)U p (z,y), 



with 



A A f x f 

i P (x,y) = -g n P (x,y) + - J o dz ° J dz i P {x,z)n f 

Kf(x, y) = -\(f 2 (x, y) - ±p 2 (x, y)) , 
Rp(x,y) = -F(x,y)p(x,y). 



{z,y), 



(13) 
(14) 

(15) 
(16) 



6 



Using (8) in (13) and (14), one observes that each term of the infinite re- 
summation of loops contained in Ip and I p contributes with the same order 
in A. To describe this regime and the late-time behavior a nonperturba- 
tive approximation has to be employed, such as the 2PI 1/iV-expansion at 
NLO. 6 ' 7 

In the fluctuation dominated regime, for t > i n onpcrt, the dynamics is 
characterized by a slow drifting of modes 6 , which suggests (cf. Ref. 17 ) the 
applicability of a first order gradient expansion with respect to the center 
coordinate X = (x + y)/2. MO.n.^.is Exploiting the effective loss 6 of the 
dependence on the initial time for sufficiently large x° and y°, we send 
t — > — oo in Eqs. (9) and (13). Using the retarded and the advanced 
propagator 

G R (x,y) = e(x°-y ) p(x,y) , 

G A (x,y) = -&(y°-x ) p(x,y) 

as well as corresponding definitions for the retarded and advanced self- 
energies E/^a and the resummation functions Ir^a, one can then also send 
the upper limits to infinity and Fourier transform with respect to the rela- 
tive coordinate s = x — y. For the two-point functions we write 

F(X,k) = J d 4 se iks F(x+^,X- S -) , 

Q{X,k) = -i J &*se iks p(x+ S -,X- S ^ . 

Here we introduced a factor — i in the definition of the Wigner transformed 
spectral function to obtain a real g (X, k). The advanced and the retarded 
propagators satisfy G* R (X, k) = Ga(X, k). Again, this property also holds 
for the corresponding self-energies. With the notation 

—jLf(X,k-q)g{X,q) 
and the definition of the Poisson brackets 



{f(X,k);g(X,k)} 



PB 



d k J(X,k)\ [dx»g(X,k) 
d x .f(X,k)] \d k ~g{X,k) 



one finds from a first-order gradient expansion of Eqs. (9) - (16) the kinetic 
equations: 



7 



2k^d x » + {d x .M 2 (X)) fc J F {X, fc) 

= P {X, fc) E e (X, fc) - S F (X, fc) q (X, fc) (17) 



[2fc^ x , + (a»M 2 (X)) d K ] q (X, fc) (18) 
= \t e (X, fc) ; Re G R (X, k)}^ + {Re S k (X, fc) ; g (X, fc)} . 
Here the Wigner transformed self-energies are given by 



S F (X,fc) = -A/(F*/ F )(X,fc) + i(^*7 e )(X,fc)j , (19) 



v , ( x. > . = - A ( (f * / e ) (x, fc) + (e * j F ) (x, fc) 



(20) 



and the equations for the resummation functions read 

i e (x,k) (g a *f) (X,k) 

- % - [l R (X, fc) ; n e (X, k)) pB + l - [l e (X, fc) ; (g a * f) {X, fc)} j 



I e (x,k) = ~[u e (x,k) i-i R (x,k) 



PB 



(21) 



7f (X, fc) = - ^ (n^ (*, fc) [l - I R (X, fc)] + I F (X, fc) (g a * P) (X, fc) 

{l R {X, k)-,u F {x, k)} pB + 1 {/ F (X, fc) ; (g a * f) (X, fc)}^ . 

(22) 

Here flp te (X,k) are the Wigner transforms of the functions defined in 
Eqs. (15) and (16). At this order of the gradient expansion the retarded 
propagator fulfills the algebraic equation 

-l 



G R (X, fc) = -k 2 + M z (X) + E R {X, fc) 



(23) 



Similarly the retarded self-energy and the retarded resummation function 
satisfy 
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t R (X, k) = — [ (F * I R ) (X, k) + (Gr * I F ) (X, k) 



(24) 



and 



I R (X, k) = -[l-I R {X, k)] (g r * f) (X, k) 



PB 



(25) 



We emphasize that Eqs. (17) - (25) represent a closed set of equations. In 
their range of applicability these equations have the advantage that they 
do not involve an integration over time history. This is important for an ef- 
ficient description of the late-time behaviour of quantum fields. A detailed 
discussion will be presented in Rcf. 8 . 

We thank Julien Serreau for collaboration on related work 2,7 and many 
discussions. 
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